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We compute A-jy^ for two flavors of light dynamical quarks using non-perturbatively 0{a) improved Wilson 
fermions. We improve on a recent calculation by employing Pade-improved two-loop and three-loop perturbation 
theory to convert the lattice numbers to the MS scheme. 



1. THE METHOD 

The running of the strong coupling constant as 
is parameterized by the A parameter. In a recent 
letter 0| we have computed A-p^ from the av- 
erage plaquette and the scale parameter tq. An 
essential step in our calculation was the conver- 
sion of the boosted coupling g^j to the renormal- 
ized coupling gj^s- This was done in two-loop 
tadpole improved perturbation theory, which re- 
sulted in an (estimated) systematic uncertainty 
of the order of the statistical error. 

We can do better now. In this talk we shall ap- 
ply Pade techniques to approximate an unknown 
three-loop coefficient in the conversion formula. 
This reduces the systematic uncertainty to an in- 
significant amount, as we shall see. 

The calculation of proceeds in three steps: 



First we compute the average plaquette P 
(1/3) {TrU\j) = Uq from which we obtain g\j, 
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where g{a) is the bare lattice coupling. Note that 
the non-perturbative contribution to the plaque- 
tte, i.e. the difference between the lattice number 
and the perturbative series, is less than lO^'^ 

In the second step we compute gjjg- from g^j. 
This is the main task. To three loops we obtain, 
following the notation of ||l| , 
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Figure 1. The four- loop [3 function for the MS 
scheme compared with various approximations. 

where — tt—pi. We know aU coefficients except 
for t^. (Sec for references.) The latter can 
be computed from the difference of the four-loop 
coefficients, b^^^ — b^, of the (3 function 

0''i9s) = -9Ubo + h,gl + big% + hig% • ■ •),(3) 
where S = MS, □ denotes the scheme. We find 
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We know b^^^ , but we do not know In the 
MS scheme the four-loop /3 function turns out 
to be very well reproduced by the Pade-improved 
three-loop P function 
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In Fig. 1 we compare (||) with the four-loop ex- 
pression. We see that in the region ff^j-g ^ 2, 
which is the region relevant to our calculation, 
the Pade-improved three-loop result differs by less 
than 0.02% from the four-loop result. We expect 
that the same is true for the boosted /3 function 
which converges at a similar rate as (3^^^ . 
Thus we may approximate by j^j . Expand- 
ing (H) gives us then 
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We can improve the convergence of the series (g) 
further by re-expressing it in terms of the tadpole 
improved coefficients 
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csw = cswUq, am — am/UQ. 
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This amounts, for example, to replacing every 
csw by 

csw (1 + ^Pi 9a + ^ (P2 - ^pf) 9n) ■ (8) 

We have already seen in 0] that tadpole im- 
provement reduces the coefficients of the two-loop 
fermionic contribution. We find that tadpole im- 
provement is equally successful in reducing the 
coefficients of the three- loop contributions. 

In the final step we compute A' 
using the renormalization group equation 
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For l3^^^ we use the Pade-improved four-loop (3 
function 
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which is a better approximation than the four- 
loop (3 function. The result should be indepen- 
dent of /Lt if our approximations of (3^^^^ in (^ 
and (||) are consistent. 

2. RESULTS 

We have varied /i between 4 and 10 GeV and 
found indeed that Aj^g- changes by a fraction of 
a percent only. Finally we have chosen a/i = 2.5, 
which corresponds to /i « 5 GeV. We may esti- 
mate the error arising from our approximation (^) 
by replacing 63-'^'^ by as well. We found 

that this changes our results by at most 0.05%. 

To convert our numbers to a physical scale we 
have used the force parameter rg . The results are 
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Table 1 

The results for A- 
simulation. 
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ro, together with rg/a, P and the quark masses am, and the parameters of the 



presented in Table 1. The numbers have increased 
by w 2% relative to our previous values [|l| , which 
is well inside the systematic error that was esti- 
mated in that paper. We fit our data by 
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^Ms''o ^A{l + Bam)[l + Cmrn) + D{a/rnf .{11) 

(For details of the fit see [0.) In Fig. 2 we show 
Ap-g ro(m = 0) = Kj^rQ — A{Bam + Cmro + 
BCammro) and Ajg^ro(a = 0) = (Ap-g-ro — 
-D(a/ro)^)/(l + Bam) which, if the fit is success- 
ful, should collapse the data points onto a single 
line each. We see that this is indeed the case. We 
also see that the largest uncertainty arises from 
the continuum extrapolation, while the data de- 
pend only weakly on the quark mass. In the chiral 
and continuum limit our fit gives 

Ajjgro = 0.553(34) (12) 

{— A). Using ro — 0.5 fm, we finally obtain 

(13) 



a5^2 = 218(13) MeV. 

This is to be compared with our previous result 
^w' = 217(16)(ll)MeV. 

3. CONCLUSIONS 

Our calculation faces two problems. One is the 
conversion from the boosted to the MS scheme. 
The other is the extrapolation to the chiral and 
continuum limit. We have solved the first prob- 
lem, and with time and more data we hope to 
improve on the extrapolation. 



Figure 2. Ajjgro{m = 0) (top) and Ajgg-ro(a 
0) (bottom) together with the fit. 
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